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We have observed exclusive 77 production in proton-antiproton collisions at = 1.96 TeV, using 
data from 1.11 ± 0.07 fb~^ integrated luminosity taken by the Run II Collider Detector at Fermilab. 
We selected events with two electromagnetic showers, each with transverse energy Et > 2.5 GeV 
and pseudorapidity |r;| < 1.0, with no other particles detected in —7.4 < rj < +7.4. The two showers 
have similar Et and azimuthal angle separation A(j> ~ vr; 34 events have two charged particle tracks, 
consistent with the QED process pp ^ p + e^e~ +p by two-photon exchange, while 43 events have 
no charged tracks. The number of these events that are exclusive ttV is consistent with zero and 
is < 15 at 95% C.L. The cross section for pp — s> p + 77 + p with |t?(7)| < 1.0 and Erij) > 2.5 GeV 
is 2.48 (Stat) l«:«(syst)pb. 

PACS numbers: 12.38.Lg, 12.40.Nn, 13.85. Qk, 14.80. Bn 



In proton- (anti)proton collisions, two direct high-i?T 
photons can be produced at leading order by qq — ^ 77 and 
by gg — > 77 through a quark loop. In the latter case it is 
possible for another gluon exchange to cancel the color 
of the fusing gluons, allowing the (anti)proton to emerge 
intact with no hadrons produced. For pp collisions, this 
is the "exclusive" process — > p + 77 + p, for which the 
leading order diagram is shown in Fig. la [1, 2]. The out- 
going (anti)proton has nearly the beam momentum, and 
transverse momentum pT ^ 1 GeV/c, having emitted a 
pair of gluons in a color singlet. There is a pseudorapid- 
ity gap At] > 6 adjacent to the (anti)proton. In Rcgge 
theory this is difFractivc scattering via pomeron [3, 4], P, 
exchange. The cross section for 177(7) | < 1-0 ^ind trans- 
verse energy Et{j) > 2.5 GeV is predicted [■'3, G] to be 
o'(77)oxciusivo ^ 0.2 - 2 pb, depending on the low- a; (un- 
integratcd) gluon density. Additional uncertainties come 
from the cross section for 3 + <?—>■ 7 + 7, the probabil- 
ity that no hadrons are produced by additional parton 
interactions (rapidity gap survival factor and Sudakov 
suppression [7]), and the probability that neither pro- 
ton dissociates (e.g., p— >p tt+tt^) [5]. The calculation 
is also imprecise because of the low Q^, the squared 4- 
momentum transfer. The total theoretical uncertainty 



on the cross section can be estimated to be a factor ^3 
[s]. Apart from its intrinsic interest for QCD, the pro- 
cess tests the theory of exclusive Higgs boson produc- 
tion [1, 2, 5, 8-13] p-\-p^p + H +p, Fig. lb, which may 
be detectable at the LHC. The leading order processes 
gg — 77 and gg^ H are calculable perturbatively, but 
the more uncertain elements of the exclusive processes 
(mainly the unintegrated gluon densities, the Sudakov 
suppression, and the gap survival probability) are com- 
mon to both (see Fig. 1). For a 120 GeV standard model 
Higgs boson the exclusive cross section at ^/s = 7 TeV is 
3 fb with a factor ^3 uncertainty [n]. 

Processes other than gg — ^ 77 can produce an exclusive 
77 final state. Contributions from qq — > 77 and 77 — > 77 
are respectively < 5% and < 1% of gg^^^ [5]. Back- 
grounds to exclusive 77 events to be considered are 7r°7r° 
and 7777, with each meson decaying to two photons, of 
which one is not detected. We also consider events where 
one or both protons dissociate, e.g., p— >p7r+7r~, to be 
background. These backgrounds are small. 

We previously published a search for exclusive 77 pro- 
duction, finding three candidate events with Exi"/) > 5 
GeV and I77I < 1.0, using data from 532 pb"^ of inte- 
grated luminosity [14]. The prediction of Ref. [")] was 
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FIG. 1. Leading order diagrams for central exclusive produc- 
tion in p{p) — p collisions: a) exclusive 77 production in p — p 
collisions; b) exclusive Higgs boson production in p — p col- 
lisions. Note the screening gluon that cancels the color flow 
from the interacting gluons. 



0.8_g 5 events. Two events had a single narrow elec- 
tromagnetic (EM) shower on each side, as expected for 
77, but no observation could be claimed. This Letter 
reports the observation of 43 events with a contamina- 
tion of < 15 7r°7r° events (at 95% C.L.), after we low- 
ered the trigger threshold on the EM showers from 4 
GeV to 2 GeV and collected data from another 1.11 fb^^ 
of integrated luminosity. We used the QED process 
p + p-^p + ^*^* + p^p + e'^e~ + p in the same data set, 
for which the cross section is well known, as a check of 
the analysis. 

The data were collected by the Collider Detector at 
Fermilab, CDF II, at the Tevatron, with pp collisions at 
Vs = 1.96 TeV. The CDF II detector is a general pur- 
pose detector described elsewhere [15]; here we give a 
brief summary of the detector components used in this 
analysis. Surrounding the beam pipe is a tracking sys- 
tem consisting of a silicon microstrip detector, a cylindri- 
cal drift chamber (COT) [16], and a solenoid providing 
a 1.4 Tesla magnetic field. The tracking system is fully 
efficient at reconstructing isolated tracks with pr > I 
GeV/c and [ryl < 1. It is surrounded by the central and 
end-plug calorimeters covering the range \r]\ < 3.6. Both 
calorimeters have separate EM and hadronic compart- 
ments. A proportional wire chamber (CES) [17], with 
orthogonal anode wires and cathode strips, is embed- 
ded in the central EM calorimeter, covering the region 
of [77) < 1.1, at a depth of six radiation lengths. It allows 
a measurement of the number and shape, in both rj and 
azimuth (f), of EM showers (clusters of wires or strips). 
The anode-wire pitch (in (j)) is 1.5 cm and the cathode- 
strip pitch varies with rj from 1.7 cm to 2.0 cm. The CES 
provides a means of distinguishing single photon show- 
ers from 7r° — >■ 77 up to i?T(7r°) ^ 8 GeV. The region 
3.6 < \ri\ < 5.2 is covered by a lead-liquid scintillator 
calorimeter called the Miniplug [!>•]. At higher pseudo- 
rapidities, 5.4 < [77] < 7.4, scintillation counters, called 
beam shower counters (BSC-1/2/3), are located on each 
side of the CDF detector. Gas Cherenkov detectors, with 
48 photomultiplicrs per side, covering 3.7 < \ri\ < 4.7, de- 



tect charged particles, and were also used to determine 
the luminosity with a 6% uncertainty [19]. 

The data were recorded using a three-level on-line 
event selection system (trigger). At the first level we 
required one EM cluster with Et > 2 GeV and \ri\ < 2.1 
and no signal above noise in the BSC-1 counters = 
5.4 — 5.9). This rapidity gap requirement rejected a large 
fraction of inelastic collisions as well as most events with 
more than one interaction (pileup). A second EM clus- 
ter with similar properties was required at level two. A 
level three trigger selected events with two calorimeter 
showers consistent with coming from electrons or pho- 
tons: i.e., passing the requirement (cut) that the ratio 
of shower energy in the hadronic (HAD) calorimeter to 
that in the EM (HAD:EM) be less than 0.125, and that 
the signal shape in the CES is consistent with a single 
shower. 

We now describe the offline selection of events, with 
two isolated EM showers and no other particles except 
the outgoing p and p, which were not detected. Two cen- 
tral, \r]\ < 1, EM showers were required with Et > 2.5 
GeV to avoid trigger threshold inefficiencies. The energy 
resolution is dE/E ^ 8% from test beam studies and in 
situ p/E matching for electrons. A refined IIAD:EM ra- 
tio cut of < 0.055 -I- 0.000 45i? was applied, as well as an 
acoplanarity cut of [tt — A(j>\ < 0.6. The trigger selection 
efficiency for single photons was measured using data col- 
lected with an interaction trigger (minimum bias). The 
BSC-1 gap trigger was taken to be 100% efficient as the 
BSC-1 trigger threshold was clearly above the noise level 
and the offline selection criteria. We measured an overall 
trigger efficiency of Etiig = 92% ± 2%(syst). A weight- 
ing process was necessary due to the different slope in 
Et of the minimum bias probe data compared to the 
signal. The trigger efficiency did not show any ry or <^ 
dependence for [77] < 1. Monte Carlo signal simulation 
data samples were generated using the SUPERCHIC pro- 
gram (version 1.3) [11, 20] based on recent developments 
of the Durham KMR model [2]. The Monte Carlo sam- 
ples were passed through a simulation of the detector, 
CDFSIM 6.1.4.m including GEANT version 3.21/14 [21]. 
The systematic error was estimated by using the bin- 
wise uncertainty of the efficiency in the weighting pro- 
cess of the signal Monte Carlo sample. Taking into ac- 
count a combined detector and offline reconstruction ef- 
ficiency of eroc = 55% ± 3%(syst), and a photon iden- 
tification efficiency of eid = 93% ± l%(syst), we ob- 
tained a photon-pair efficiency £pho = ^trig * ^rcc * ^fd — 
40% ± 3%(syst). The systematic uncertainties of the re- 
construction and identification efficiency were estimated 
by shifting kinematical input parameters over a reason- 
able interval motivated by the dominating EM-energy- 
scale uncertainty [22]. The offline selection then required 
that no activity other than these two showers (or clusters 
of showers) occurred in the entire detector, ]?7] < 7.4. We 
used the same procedure as in our earlier study of exclu- 
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events [23], searching all the calorimeters for 
any signal above noise levels, determined using noninter- 
action events; 99.2% of such events have no tower (out 
of 480) with Et > 125 MeV. We also required the CLC 
counters and the more forward BSC counters to have 
signals consistent with only noise. Events triggered only 
on a bunch crossing (zero bias) showed that the exclu- 
sive efficiency, Sexci, defined as the factor to be applied 
to the delivered luminosity to account for the require- 
ment of no pileup, is Sexci = 6.8% ± 0.4%(syst). The 
probability P(0) of a zero-bias event satisfying all the 
exclusivity cuts, i.e., having no detected inelastic inter- 
action, is P(0) = j4exp(— n) = Aexp(— I/xCvis),where 
is the single bunch crossing luminosity (cm~^) and dvis 
is the visible cross section; dvis = find if cveryinelastic 
collision is detected. We find (Tvis = 67 ± 6 mb. In the 
absence of noise (above our chosen thresholds) A = 1.0; 
we find A = 0.98 ±0.02. We checked that the rate of can- 
didate events, corrected for the exclusive efficiency, was 
constant during data taking (one year). The systematic 
uncertainty was estimated using the spread in slope pa- 
rameters from fits to data in different time periods. 



TABLE I. Summary of parameters used for the measurement 
of the exclusive photon-pair cross section for Et{j) > 2.5 GeV 
and |?7(7)| < 1.0. Values for the e"^e~ control study are also 
given. Note that b/g stands for background. 



Integrated luminosity C;„t 
Exclusive efficiency 



1.11 ±0.07 fb"" 
0.068 ± 0.004 (syst) 



Exclusive 77 



Events 

Photon-pair efficiency 
Probability of no conversions 

b/g (events) 
Dissociation b/g (events) 



43 

0.40 ± 0.02 (stat) ± 0.03 (syst) 
0.57 ± 0.06 (syst) 
0.0, < 15 (95% C.L.) 
0.14 ± 0.14 (syst) 



Events 

Electron-pair efficiency 
Probability of no radiation 
Dissociation b/g (events) 



34 

0.33 ± 0.01 (stat) ± 0.02 (syst) 
0.42 ± 0.08 (syst) 
3.8 ± 0.4 (stat) ± 0.9 (syst) 



The selection of 81 events passing all cuts was made 
without reference to the track detectors. We found that 
34 have exactly two oppositely charged tracks, 43 have 
no tracks in the COT, and four are in neither class. Vi- 
sual inspection of the latter showed that two had photon 
conversions, and two were likely to be e"''e~ events with 
bremsstrahlung. These numbers are consistent with ex- 
pectations from the detector simulation. The tracks in 
the 34 two-track events agree in all aspects with the QED 
process p+p p+e'^e^ +p via two virtual photons, pre- 
viously observed in CDF [23, 24]. The calorimeter shower 
energies are consistent with the momenta measured from 
the tracks. Kinematic distributions, after detector sim- 
ulation, are as expected. The mass M{e^e~) distribu- 
tion is presented in Fig. 2a, together with the QED pre- 
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FIG. 2. The e^e~ candidates: invariant mass distribution 

(a) . The two-photon candidates: invariant mass distribution 

(b) , \tt ~ A<^| distribution (c), and pr distribution of the two 
photons (d). All error bars are statistical. The MC predic- 
tions for 77 are normalized to data. The QED prediction for 
e^e^ is normalized to the delivered luminosity and efficien- 
cies. The MC samples for the QED process were generated 
with the LPAIR program [25]. 



diction normalized to the delivered luminosity and effi- 
ciencies, showing that the cross section agrees with the 
QED prediction in both magnitude and shape. We mea- 
sured a cross section of CTg+e-, exciusive(|^(e)| < 1, Erie) > 
2.5 GeV) = 2.88 toUistat) ± 0.63(syst) pb, compared to 
3.25±0.07pb (QED, [2.j]). The systematic uncertainties 
for the QED study are mostly identical to the photon 
case. Distinct from photons, electrons leave tracks in the 
tracking detectors and may radiate. The systematic un- 
certainty on the radiation probability was estimated by 
varying the exclusivity cuts by ±10%. This e+e~ sample 
provides a valuable check of the exclusive 77 analysis. 

The 43 events with no tracks have the kinematic prop- 
erties expected for exclusive 77 production [20]. In par- 
ticular the M(77) distribution [Fig. 2b] extending up to 
15 GeV/c^ is as expected, as well as the acoplanarity 
TT — A(/)(77) [Fig. 2c] and the 2- vector sum of [Fig. 2d]; 
in these plots [unlike Fig. 2a] the SUPERCHIC Monte Carlo 
prediction is normalized to the same number of events as 
the data. An important issue is whether some of these 
events could be 7r°7r°, rather than 77. Note that 77r° 
events are forbidden by C parity. The CES chambers 
give information on the number of EM showers. The 
minimum opening angle AOmin between the two photons 
from 7r° decay is 2tan~-'^ (^^) ~ P^^^ ~ ^ 
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GeV, well separated in the CES chambers, which have a 
granularity < 0.5°. A tt" can fake a 7 only if one photon 
ranges out before the CES, or falls in an inactive region 
(8%) of the detector. All of the 68 e± events m our sam- 
ple, with similar energies, had matching showers in the 
CES chambers. A GEANT [21] simulation predicts the 
probability that a photon in our energy range produces 
a shower to be > 98.3%. We summed the number of re- 
constructed CES showers in the event, mostly 2 or 3 as 
shown in Fig. 3 (left). The distribution agrees very well 
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FIG. 3. Estimate of tt tt background fraction in the candi- 
date sample. Distribution of reconstructed CES showers per 
event for data compared to 77 and n'^n" Monte Carlo simu- 
lations (a). Background fraction estimate using Pearson's 
test to fit the composition hypothesis to the data distribution 
(b). 



with the 77 simulation, and strongly disagrees with the 
TT^TT^ simulation. Fitting to the sum of the two compo- 
nents gives a best fit to the fraction F(7r°7r°) = 0.0, with 
a 95% C.L. upper limit of 15 events. Since obtaining 
this result, a new calculation of exclusive tt'^tt'^ produc- 
tion [26] predicts aexdiT^^T^") 6 - 24 fb for EriTr") > 2.5 
GeV and < 1.0, < 0.01 of our measured exclusive 77 
cross section. In the cross section calculation we take 
this background to be zero. Exclusive 7777 production is 
also expected to be negligible. The only other signif- 
icant background could be undetected proton dissocia- 
tion, about 10% for the QED e+e~ process but <1% 
for P+P— >■ 7 -|- 7 [5, 27, 28]. The cross section for both 
photons with Et{j) > 2.5 GeV and [77(7)1 < 1-0 ^-nd no 
other produced particles is given by: 



a. 



77, exclusive 



N (candidates) — N (background) 

H'in^.S Qxcl 



where e is the product of the trigger, reconstruction, 
identification, and conversion efficiencies (22.8%) in Ta- 
ble I. The systematic uncertainty on the conversion prob- 
ability was estimated by varying the exclusivity cuts by 
±10%. We find (T77,exei (|r?(7)| < l,ET(-f) > 2.5 GeV) = 
2.48io;35(stat)lo;5i(syst) pb. The theoretical predic- 
tion [11] is strongly dependent on the low- 2; gluon den- 
sity, having central values 1.42 pb (mstwOSlo) or 0.35 



pb (mrst99), with other uncertainties estimated to be 
a factor of about ^3 [28]. A comparison of our mea- 
surement with the only theoretical prediction available 
to date is shown in Fig. 4. The rates of e+e~ and 77 
events with ET(e/j) > 5 GeV are consistent with those 
in our earlier studies [14, 23]. 
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FIG. 4. Comparison of the measured cross section for the 
exclusive 77 production in pp collisions at ^/s = 1.96 TeV 
with theoretical predictions [I I]. 



In conclusion, we have observed the exclusive produc- 
tion of two high-Ex photons in proton-antiproton colli- 
sions, which constitutes the first observation of this pro- 
cess in hadron-hadron collisions. The cross section is in 
agreement with the only theoretical prediction, based on 
g+g — >■ 7+7, with another gluon exchanged to cancel the 
color and with the p and p emerging intact. If a Higgs bo- 
son exists, it should be produced by the same mechanism 
(see Fig. 1), and the cross sections are related. 
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